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Abstract 

A survey is given of the present state of the art in studying Bose-Einstein condensation 
of dilute atomic gases. The bulk of attention is focused on the principal theoretical prob- 
lems, though the related experiments are also mentioned. Both uniform and nonuniform 
trapped gases are considered. Existing theoretical contradictions are critically analysed. A 
correct understanding of the principal theoretical problems is necessary for gaining a more 
penetrating insight into experiments with trapped atoms and for their proper interpretation. 
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1 Introduction 



Bose- Einstein condensation of dilute atomic gases is intensively developing field of research. There 
exist several reviews and books on the subject, for instance [1-5]. However, both experiment and 
theory in this field develop so fast that voluminous works are not able to follow the front line of 
research. Moreover, the amount of information in the field is so huge that even detailed reviews do 
not include all aspects of the problem. The aim of this survey is to present a brief description of 
the modern state of the art related to Bose-Einstein condensation of uniform and trapped dilute 
gases. The main concern will be given to the principal problems in Bose-Einstein condensation. 
Some of these problems gave rise to controversies. It is, therefore, important to present a critical 
analysis of the basic theoretical foundations, whose correct understanding is compulsory for the 
proper interpretation of experiments with trapped Bose-condensed atoms. 

Nowadays, there are several tens of laboratories in different countries, where Bose-Einstein 
condensates are routinely produced. Initially, the condensation of atoms with nonzero spins has 
been realized, with the total atomic spin F = S -|- L -|- I corresponding to a boson, where S is 
electron spin; L, electron angular momentum; and with I denoting the nuclear spin. The following 
dilute atomic gases have been condensed: ^^Rb [6], ^Li [7], ^^Na [8], ^^Rb [9], ^^K [10], ^^^Cs [11], H 
[12], and metastable triplet ''He [13,14]. Recently, the spinless atoms of ^''^Yb [15] were condensed. 
There are indications of producing ultracold condensed molecules ^^Rb2 [16] formed, by means 
of Feshbach resonance, in the Bose-Einstein condensate of ^^Rb. The condensation of molecules 
^Li2 [17], from the degenerate Fermi gas of ^Li has been reported. More fully the condensation of 
molecules will be detailed in the last section below. 

In this review, only the systems with Bose-Einstein condensates will be considered. There 
exists now an intensive activity of studying degenerate trapped Fermi gases. These, however, 
will not be touched here, except the case when Fermi atoms form Bose molecules that undergo 
condensation [17]. In the frame of a brief survey, it is impossible to cover all topics related to cold 
trapped atoms. So, the emphasis will be given to the principal theoretical points in Bose-Einstein 
condensation of dilute atomic gases. 



2 Condensation Temperature 

The infiuence of atomic interactions on the condensation temperature is a rather nontrivial prob- 
lem, because of which much efforts have been devoted to finding the dependence of the conden- 
sation temperature Tc on effective interaction parameters. The interaction in dilute gases is well 
characterized by the Fermi contact potential 

$(r) = 47r ^ S(r) , 
m 

in which is the s-wave scattering length and m, atomic mass. The gas is dilute when \pa^\ <S 1, 
where p = N/V is the number density. The problem to be solved is how Tc depends on pa^. 

This problem turned out to be complicated even for homogeneous dilute gases, since Bose- 
Einstein condensation is a second-order phase transition, which is governed by long-range correla- 
tions. Because of these correlations, simple mean-field approximations are not always applicable. 
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In an ideal homogeneous Bose gas, the condensation temperature is 



^0 = 



mk 



B 



am 



2/3 

(1) 



where ({■) is the Riemann zeta function, with C(3/2) = 2.612. What happens when in the ideal gas 
the atomic interactions arc being switched on, while the density is kept fixed? The corresponding 
change of the condensation temperature can be characterized by the relative difference 

at; i; 

= — - 1 . (2) 

To To ^ ^ 

This quantity has been calculated by a variety of methods. Historical references can be found in 
[5,18,19]. For repulsive interactions, when the scattering length is positive and changes from zero 
to a finite value a^, the generally accepted functional dependence of the relative shift (2) on the 
diluteness parameter or gas parameter 



6 = p ' tts^ — 



is given by the asymptotic expansion 

AT. 



~Ci5 + (c'2ln5 + c'2') , (3) 



as (5 — >^ 0. The first-order coefficient ci depends on nonperturbative correlation effects. Its most 
accurate value was found by numerical Monte Carlo calculations for three-dimensional lattice 0(2) 
field theory, which gave ci = 1.29 ± 0.05 [20,21] and ci = 1.32 ± 0.02 [22,23]. The coefficient 4 
can be computed exactly using perturbation theory [24], which yields 

647rC(l/2) 

3[C(3/2)]^/3 - '^-^''^ ■ 

The remaining second-order coefficient depends on nonperturbative physics and can also be 
extracted from numerical lattice simulations of three-dimensional 0(2) scalar field theory [23,24], 
resulting in = 75.7 ± 0.4. Summarizing these results, one has 

ci = 1.32 ± 0.02 , 4 = 19.7518 , 4 = 75.7 ± 0.4 . 

There have also been many attempts to calculate the nonperturbative coefficients Ci and C2 in 
different approximations. Thus, within the frameworks of Ursell operators and of Green functions 
[18], the results are ci = 1.9 ,03 = 5.589, and C2 = 4.9, which arc rather far from the correct 
values. Renormalization-group calculations give Ci = 1.23 [25] and Ci = 1.15 [26], but C2 was not 
calculated. The 1/n-expansion in the n-component field theory has also been used, yielding in 
the first and second orders, respectively, ci = 2.33 [27] and ci = 1.71 [28], but, again, C2 was not 
computed. We shall not enumerate here many other attempts, trying to find the coefficients in 
the shift (3), but whose results arc not sufficiently reliable (see discussion in [18,29]). 

The most accurate results, obtained by means of an analytical approach, and being closest 
to the Monte Carlo numerical simulations, are due to the usage of the optimized perturbation 
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theory. The main idea of this theory, whose general formulation was first given in [30] (see further 
references in survey [31]), is the introduction of control functions governing the convergence of the 
sequence of optimized approximants. Different variants of the optimized perturbation theory are 
now widely used for various problems [31], in particular, in quantum field theory [32-38]. There 
are several ways of introducing control functions. 

A straightforward method is the introduction of control functions at the initial stage of a calcu- 
lational procedure, when they can be incorporated into the trial Hamiltonian or Lagrangian, or into 
the trial wave function or Green function, which the following perturbation theory will be based 
on [30] . Control functions are defined at each step of the employed perturbation theory so that to 
make convergent the sequence of the resulting optimized approximants. This method of introduc- 
ing control functions into the initial Lagrangian of field theory has been used for calculating the 
critical temperature shift in [39-42], where at the second, third, and fourth orders of the optimized 
perturbation theory it was found, respectively, Ci = 3.06, 2.45, 1.48 and C2 = 101.4, 98.2, 82.9. 
These sequences seemingly converge to the correct numerical Monte Carlo values. The problem of 
convergence was specially studied in [43-46], where the optimized perturbation theory was com- 
bined with 1/n-expansion. The series are shown to converge exponentially. The final estimates 
for the coefficients are Ci = 1.19 and d!^ = 84.9. Thus, the optimized perturbation theory gives 
the results for both coefficients ci and C2 that exponentially converge to the Monte Carlo values 
[39-46]. 

It is possible to notice that the introduction of control functions into the initial Lagrangian or 
Hamiltonian leads to a special scaling of the calculated quantities. Using this, a way of introducing 
control functions, in the frame of the optimized perturbation theory, can be done at the final stage, 
when one, first, derives divergent series and then transforms them, with control functions being 
incorporated in the process of this transformation. A particular transformation method has been 
suggested by Kleinert [47]. This method was invoked for calculating the coefficient ci [48-50]. 
The calculation, first, starts with deriving loop expansions for (3-|-l)-dimensional field theory. 
The highest expansion, derived for ci, is the seventh- loop expansion [50] 

ci{g) ~ aig + a2g^ 03/ 04/ a-^g^ 

in powers of an effective coupling parameter, assumed here to be asymptotically small, g ^ The 
actual ci corresponds to the limit ci(oo) g ^ 00. As is clear, the above small-^f expansion has 
no sense for g ^ 00. To find an effective limit ci(cxo), one has to renormalize the above expansion. 
Such a renormalization can be done by introducing control functions by means of the Kleinert 
transformation [47] and by defining these control functions from an optimization procedure. This 
way also leads to an exponentially fast convergence of the sequence of optimized approximants 
[47]. The coefficients ak depend on the number of components n. For the two-component field, 
corresponding to a Bose system, one has [50] 

oi = 0.223286 , 02 = -0.0661032 , 03 = 0.026446 , 

a4 = -0.0129177 , 05 = 0.00729073 (n = 2) . 

Employing the Kleinert method, it is necessary to assume that Ci(oo) is finite. This limit also 

essentially depends on the choice of the Wegner exponent, characterizing the approach to the 
limit Ci(cxo). Kastcning [49,50] calculated Ci = Ci(oo) for two different values of the Wegner 
exponent. When the latter is defined self- consistently, from the same expansion, then ci = 1.376; 
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but if this exponent is fixed, being chosen from some additional arguments, then ci = 1.611. The 
average value C\ = 1.269 seems to be close to the Monte Carlo results. To check the method, 
Kastcning [50] has also calculated Ci(cxd) for one-component and four-component fields, for which 
Monte Carlo simulations are available as well [51]: Thus, for n = 1, Ci = 1.09 ± 0.09 and for 
n — A, ci — 1.59 ± 1.10. For n — 1, the expansion coefficients are 

ai = 0.334931 , = -0.178478 , = 0.129786 , 

04 = -0.115999 , og = 0.120433 (n = 1) . 

The seven-loop optimized approximant [50] is ci = 1.171 for the self-consistent choice of the 
Wegner exponent and ci = 0.973, if the latter is fixed. The average ci = 1.072 is close to the 
Monte Carlo value. In the case of the four-component field, 

oi = 0.167465 , 02 = -0.0297465 , 03 = 0.00700448 , 

04 = -0.00198926 , 05 = 0.000647007 (n = 4) . 

The two choices of the Wegner exponent lead [50] to Ci = 1.648 and ci = 1.435, whose average 
Ci = 1.542 is again in agreement with the Monte Carlo result. The coefficient C2 was not calculated 
by this method. 

The accuracy of the optimized perturbation theory can be essentially improved by introducing 
several control functions. This has recently been done [52] for calculating the coefficients in the 
temperature shift (3). For the coefficient ci, the seven- loop expansion was used, yielding ci = 1.3 
and for C2, the six-loop expansion was employed, giving C2 = 73.46. 

The results of all reliable calculations show that the critical temperature shift (3) is positive and 
increases with the increasing diluteness parameter. That is, the temperature of Bose condensation 
grows under increasing density or under stronger interaction. Such a behaviour may look rather 
strange, though one should not forget that this happens in the dilute limit, when the diluteness 
parameter 5 = p^^^ds is small. When this parameter is large, as for liquid helium, Monte Carlo 
simulations show that the critical temperature decreases with growing density [53]. This is in 
agreement with the existence of the Bogolubov depletion, according to which the condensate 
density po zero temperature is given by the equation 

where 5 = p^^^Ug is the diluteness parameter and 

97r\^/3 
5^={ — ] = 0.761618 . 
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The Bogolubov depletion formula (4) demonstrates that there exists a critical value Sc of the 
diluteness parameter 6, such that the condensate density po becomes zero at the critical tempera- 
ture Tc = 0. The Bogolubov approximation, as is known, is valid for weakly nonideal gases, when 
5 <^ 1. Therefore the critical diluteness parameter 5c, when the complete depletion occurs, may 
not be very accurately estimated in this approximation. However Monte Carlo simulations for a 
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wide diapason of S show that the Bogolubov depletion formula is qualitatively correct for 5 ~ 1 
as well, though it is quantitatively accurate up to 5 0.1 [54]. 

In this way, the Bosc-Einstcin condensation temperature Tc as a function of the diluteness 
parameter 6 behaves as follows. First, when 6 is asymptotically small, grows with increasing S 
according to Eq. (3). But then it turns down, becoming a decreasing function of 6. And goes to 
zero at 5 ~ 1. This type of reentrant behaviour in the whole range of 5 has been studied [55] in the 
Bogolubov approximation. The maximum T^^^ of Tc was found to occur at S — Smax ~ 0.25 Sc, 
so that 

fa 1.1 , Smax ~ 0.2 . (5) 

-to 

Extracting from [54] the behaviour of as S tends to 6c from below, we get 



-f ~ 2.582^5, -S (S^Sc). (6) 
-to 

Recall that Eqs. (3) and (6) have to do with a homogeneous three-dimensional system. 

For trapped gases, the situation is quite different. Let the trapping be realized by the harmonic 
potential 

Strictly speaking, there is no sharp phase transition for a finite number of atoms N. To rigorously 
define the Bose-Einstein condensation, one should consider N ^ oo. For the gas of Bose atoms, 
trapped in the harmonic potential, the thermodynamic limit is defined [2,4] as 

A" — > oo , a;o ^ , A^a;Q — > const , 

which follows from the definition of thermodynamic stability [2] , and where the effective frequency 
is the geometric average 

In this limit, the Bose-Einstein condensation temperature of an ideal gas is 

1/3 

(7) 



1 n 



N 



C(3) 



For other types of the trapping potential, the condensation temperature will be different [2]. 
Taking into account the final-size corrections, one has [4] the ideal-gas condensation temperature 



where uj is the arithmetic average 



, _ C{2)Tiuj 
203)^' 



- 1 / N 



Similarly to the homogeneous case, it is possible to introduce the diluteness parameter for 
trapped atoms, whose definition, however, is different. For the harmonically trapped gas, the 
diluteness parameter is given by the ratio 

6 = as/Xo (9) 
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of the scattering length ag to the de Boghe wavelength at the transition temperature 



\ I. rr = V 27r lo 



N 



with Iq = ^^Ti/mujQ being the oscillator length. Switching on interatomic interactions [56] changes 
the condensation temperature to 



C(2)o7 



To 2C(3)cuo 



C(3) 
N 



1/3 



+ Ci5 + (c'2ln5 + 4V^ (10) 



where 5 is asymptotically small. In the case of a uniform Bose gas, the first-order shift in T,., related 
to the coefficient Ci, is sensitive to critical fiuctuations and so is nonperturbative. In contrast, the 
first-order shift for a gas trapped in a harmonic potential is calculable using perturbation theory 
[4]. The second-order logarithmic coefficient c'a is also calculable in perturbation theory, but the 
coefficient 03 is not [56]. The coefficient 03 can be calculated by involving Monte Carlo lattice 
simulations in three-dimensional 0(2) field theory. The results for these coefficients [56] are 

ci = -3.426032 , 4 = - = -45.856623 , 4' = -155. 

The very first thing that catches the eye, comparing the critical temperature shifts (3) and (10), 
is that for a homogeneous gas all coefficients Ci, c^-, and c'2 are positive, while for a trapped 
gas all these coefficients are negative. This makes the principal difference between the uniform 
and trapped gases. Switching on repulsive interactions in a uniform gas shifts the condensation 
temperature up, while in a trapped gas this shifts the condensation temperature down. Monte 
Carlo investigations [54], of a trapped gas at zero temperature in a wide range of the diluteness 
parameter 5 show that the condensation fraction is well described by the Bogolubov depletion 
formula (4) up to 5 ~ 0.1, but for larger 5 ~ 1, the depletion formula (4) is only qualitatively 
correct. It looks like the depletion formula (4) can be extended for high 5 as well, but with the 
critical diluteness 5c being different from that given by the Bogolubov approximation. 

As is seen, there is no the reentrant effect in trapped gases, in contrast to uniform gases, where 
the condensation temperature first increases and then diminishes as a function of the diluteness 
parameter. For trapped gases, Tc is monotonically diminishing function of the gas parameter 5. 

The formula (10) for the critical temperature of a dilute trapped Bose gas is rehable. There were 
some claims that the linear in 6 shift could be compensated by the finite-size corrections. However, 
as follows from (10), all shifts have the same negative sign, so there can be no compensation 
between different terms. In addition, formula (10) has recently been checked for ^''Rb and found 
to be in good quantitative agreement with experimental data [57]. The principal difference in the 
behaviour of the critical temperature for uniform and trapped gases can be explained by different 
properties of fluctuations in these systems. 

As will be explained in the following section, the uniform ideal gas is an unstable system, 
with anomalously strong particle fluctuations. Switching on atomic interactions stabilizes the 
system, making the fluctuations normal. As a result of such a stabilization, the condensation 
temperature, flrst, increases according to Eq. (3). But further strengthening interactions leads 
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to the condensate depletion. This is why the condensation temperature Tc as a function of the 
diluteness parameter p^^^Cg possesses a specific reentrant behaviour [55]. 

Contrary to the uniform case, the trapped ideal Bosc gas is stable, having normal fluctuations. 
Atomic interactions play only the destructive role, depleting the condensate. Therefore the critical 
temperature (10) is a monotonely decreasing function of the interaction strength characterized by 
the parameter Og/Ao defined in Eq. (9). 



3 Condensate Fluctuations 

Studying fluctuations in statistical systems is important from several points of view. For instance, 
fluctuations define the stability of the system and its way of reaching the state of thermodynamic 
equilibrium [58,59]. 

In the case of Bose-condensed systems, one usually studies the number-of-particle fluctuations 
because of their nontrivial properties. Since in literature, there are controversies related to this 
problem, it is necessary to give a careful analysis of the situation. 

The number-of-particle fluctuations are characterized by the dispersion 

A2(7V) = < > - < TV >2 (11) 
for the number-of-particle operator 

TV = y V'^(r)V'(r) dr , 

with ■0(r) being a field operator. The average number of particles A'" is a statistical average 

N ^ < N > = j p(r) dr , 

where p(r) =< ^^(r)^(r) > is the density of particles. Here we keep in mind an equilibrium state, 
because of which p(r) does not depend on time. But p(r) depends on the spatial variable r, when 
the system is nonuniform. 

The dispersion (11) is directly related to the isothermal compressibility 

"-^-fO^^Ks^L '''' 

by the equality 

whose derivation can be found, e.g., in [60]. A necessary condition for a system to be stable is the 
semi-posit iveness and finiteness of the compressibility, that is, < < oo. If the compressibility 
(12) were infinite, this would mean that an infinitesimal fiuctuation of pressure P would lead to 
an immediate collapse or explosion of the system. Therefore, in the thermodynamic limit, the 
dispersion (11) should behave as 

(TV) -const- AT {N ^ oo) . (14) 
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If A^(iV) would be proportional to A?^ in a power higher than one, this would mean that the system 
is unstable. When the stability condition (14) is satisfied, the number-of-particle fluctuations are 
called normal, but when Eq. (14) is not valid, so that the compressibility (13) diverges in the 
thermodynamic limit, the fluctuations are termed anomalous. 
The dispersion (11) can be presented as 

A' (TV) = y [R{r, r') - p(r)p(r')] dr dv' (15) 
through the density-density correlation function 

i?(r,r') = < V'^(r)V'(r)V'^(r')V'(r') > . 
The latter is related by the equation 

i?(r, = p(r)(5(r - r') + p(r)p(r')^(r, r') 
with the pair correlation function 

^ < i;\v)i,\v')i^{v')^{v) > 
p{v)p{v') 

Therefore, another representation of the dispersion (15) is 

A2(Ar) = AT + j p(r)p(r') [g{Y, r') - 1] dv dr' . (16) 

Note that the pair correlation function possesses the property g{r,r') — g{r',r). One more form 
for the dispersion (11) can be obtained by recollecting the definition of the structural factor 

S{k) ^ i I [R{r, r') - p(r)p(r')] e-^-^-^') dr dr' , 

which can also be written as 

5(k) = 1 + ^ / p(r)p(r'Mr, r') - 1] e'^-^-^') dr dr' . 

Then, from Eqs. (15) or (16), it follows that 

A^{N) = NS{0) . (17) 

All equations above are valid for any nonuniform system. In the particular case of a uniform 
system, when p(r) = p and ^'(r, r') = g{r — r'), one has the well known expression 

5(k) = l + p|[^(r)-l] e-^^-^ dr. 

Then the dispersion (16) reduces to 

A\N) = N^^l+pJ[g{r)-l]dr'^. (18) 
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We see that particle fluctuations are closely related to the behaviour of the pair correlation 
function as well as to the structural factor, which can be measured by means of light or neu- 
tron scattering [61,62]. One more useful relation comes from the expression of the isothermal 
compressibility (12) through the isothermal sound velocity s, which reads 

1 2_ 1 (dP\ 

pms^ m \op 



TN 



Then, for the zero- vector structural factor, one has 



S{0) = pkeTKr 



So that for the dispersion (17), we find 



A2(iV) = M!Ar. (19) 



Note that the same form (19) can be obtained directly from the Bogolubov theory which provides 
a good first approximation for dilute gases. Really, it is easy to show [2] that in the Bogolubov 
approximation 



.2 



-5(0 = , c = W — - pas , 

mc^ V m'^ 

with c being the sound velocity of elementary excitations. Hence, we obtain the same form (19), 
with s — c. 

Analyzing the relation of the dispersion (18) with the pair correlation function, we may recall 
that the general behaviour of the latter at large r = |r| is 

C e"""/^ 

where C is a constant, ^ is coherence length, d is space dimensionality, and r] is the index of 
anomalous dimension, < < 1. Therefore, for any finite ^, one has 5'(0) ~ const and A^TV ~ N, 
that is, the fluctuations are normal, as it should be for any stable equilibrium system. Fluctuations 
can become anomalous only when the coherence length diverges, ^ = oo, which occurs at the point 
of a phase transition. Then 5'(0) ~ iV^^"'')/^ and 

A' (TV) ~ Ar(5-^)/3 (^ = oo) . 

With the largest value 77 = 1, we have A2(iV) - N^/^. This means that the fluctuations are 
anomalous, the compressibility kt ^ 00, the sound velocity s — > 0, and the structural factor 
5'(0) — > 00. But all this is not a surprise, since the point of a phase transition is the point of 
instability. Above as well as below this point, the coherence length is finite and fiuctuations are 
normal, which tells that the system is stable. For instance, in the interacting Bose systems [63,64], 
the coherence length is ^ ~ h/ms. 

Since particle fluctuations are directly related to the behaviour of the pair correlation function, 
it is necessary to keep in mind the general properties of the latter, which follow from those of the 
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second-order density matrix [65] . The pair correlation function satisfies the asymptotic boundary 
condition 

^(r,r')^l (|r-r'|^oo) 
and the exact normahzation condition 

/^^(-')--'-. 

which for the uniform case reduces to 

If 1 
— / q(r) dr — 1 — — . 
V J ' N 

In this normahzation condition, the thermodynamic hmit is assumed. 

The fictitious appearance of anomalous fiuctuations is often due to the fact that in the process 
of calculations the general properties of the pair correlation function have been occasionally broken. 
As an illustrative example, let us adduce the following. From Eq. (18), it is clear that the major 
contribution to the integral comes from the behaviour of g{Y) at large r oo. At large distance, 
the interparticle interactions vanish. Then, using the Wick decoupling for < ip'^ {y)iP'^ {fS)ip{fd)ip{Y) , 
one could get 

< V'^(r)V^(r) >< V'HO)^(O) > + < V'Hr)V'(O) >< ip\'^)i^{v) > . 

The density of a Bose-Einstein condensate, po, can be defined as the limit | < V''''(r)'0(O) > | — > po 
for r oo. Consequently, for g(r) at large r — >■ oo, one would get 1 + N^/N"^, which would result 
in the anomalous fiuctuations A^(A^) ~ A^^q. However, the obtained limit for g{r) is in agreement 
with the asymptotic boundary condition g{r) ^ 1, as r —> oo, only for Nq = 0. Hence, actually, 
there are no anomalous fluctuations in this case. The contradiction has arisen because the Wick 
decomposition cannot be used for correlated systems with long-range order. 

The particle fiuctuations below Tq for the ideal uniform Bose gas have been studied in a 
systematic way [66,67]. In the grand canonical ensemble, the fluctuations of the occupation 
number A^^ for the operator of the numbcr-of-particles with momentum k, such that |k| > 0, 
behave as A^(iV/j) ~ N^^'^. Hence, these fluctuations are anomalous. Even more anomalous are 
the fiuctuations of the total number of particles A^(iV) = A^(Ao) = A'g, which are caused by the 
condensate fiuctuations. However, as is noted above, such fiuctuations are in agreement with the 
correct properties of the pair correlation function only when A'o = 0. Calculating the dispersion for 
the number of particles out of the condensate, Nout = J2k>o^k, one has [66,67] A\Nout) ~ A^^/s. 
In the canonical ensemble, when the total number of particles is rigorously fixed, taking into 
account the equality No — N — Nout, one gets A^(Ao) = A^(iVo„t), which is given by 

A2(7Vo)= 4.372 N''\ (20) 

where Tq is the condensation temperature (1) for the uniform ideal Bose gas. Formula (20) was 
first derived by Hauge [68]. At all finite temperatures, the condensate fluctuations are anomalous. 

Thus, the ideal uniform Bose gas is a rather pathological object, with anomalous fluctuations. 
It is easy to find its pressure [69], which for T < Tq is 
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C(5/2) 



Xt = 



\ mkBT 



(21) 



As is evident, OPq/Dp — 0, so that the comprcssibihty (12) diverges, Kt — * oo, which is in 
agreement with the anomalous behaviour of the dispersion A^(iV). Therefore, the ideal Bose 
gas is unstable below Tq, though it is stable above this temperature. For a stable equilibrium 
system, the compressibility must be finite, fluctuations normal, and also, the observable quantities, 
calculated in different statistical ensembles, must coincide in the thermodynamic limit. Hence, 
the ideal uniform Bose gas cannot exist as a stable equilibrium system below the temperature of 
condensation. 

The theory of weakly nonideal uniform Bose gas was developed by Bogolubov [70-72]. The 
dispersion A^(iV), which he calculated [71], under the conditions T < Tq and N — Nq <^ N, takes 
the form 

A'{N) = ^ iV , (22) 



where c = \JAtiTi^ pag/rn? is the Bogolubov sound velocity. This is in agreement with the general 
expression (19). The Bogolubov theory [70-72] is based on three points. First, the condensate is 
described by introducing the condensate field ipQ by means of the so-called Bogolubov prescription 
■0 = ■0o+'0 with ■00 treated as a nonoperator term and %l) = Sfc^o o>k^^^ ^ being the field operator 
of noncondcnscd atoms. Then, owing to the assumption N—Nq <C N, all terms of the Hamiltonian 
containing the factors of ip, or a^, of orders higher than two, are omitted. And finally, the remaining 
approximate Hamiltonian is diagonalizcd by means of the Bogohibov canonical transformation. 
Since only the terms of up to the second order in ip, or a^, are retained in this approximation, the 
same has to be done in calculating all observable quantities. Considering the expressions that are 
higher than quadratic in ip would go beyond the accuracy of the Bogohibov approximation and 
so would be not self-consistent. This was taken into account by Bogolubov [71] when calculating 
the dispersion (22) using the relation (18). 

The weakly nonideal Bose gas possesses the pressure [69] 



2Trh a. 



m 



^ am 



(23) 



where Pq is the pressure (21) of the ideal gas. From here, the compressibility (12) is kt — l/pmc. 
Therefore, due to Eq. (13), the particle dispersion coincides with the Bogolubov form (22). 

In this way, including atomic interactions makes the fluctuations normal, thus, stabilizing the 
system. Returning to the ideal gas would mean setting a,^ 0, because of which c and 
Kt oo, so that the system would become unstable. Equation (22), as is clear, corresponds 
to the grand canonical ensemble. In canonical ensemble, the fluctuations of the total number of 
particles are not defined, but one can find the dispersion A^(7Vo), for the number of condensed 
particles, assumed to be equal to A'^{Nout)- For a weakly excited interacting gas, condensate 
fluctuations were found to be normal, that is, A^(iVo) ~ N, both in the canonical as well as in 
the microcanonical ensembles [73]. 

In the trapped ideal Bose gas, particle fluctuations are normal, contrary to those in the uniform 
case. Politzer [74] calculated A^(7Vo) in both canonical and grand canonical ensembles for atoms 
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trapped inside a harmonically confining potential, the difference between the ensembles vanishing 
in the thermodynamic limit. The Politzer result is 

where Tq is the condensation temperature (7) for the harmonically confined ideal gas. Micro- 
canonical ensemble also displays normal condensate fluctuations [75]. Comparing Eq. (24) with 
Eq. (20), wc see that the anomalous fluctuations of the ideal uniform gas are suppressed by the 
confining potential. 

The influence of interactions in the dilute trapped gas on the size of fluctuations can be under- 
stood by remembering that these interactions diminish the condensation temperature according 
to Eq. (10). Therefore, the dispersion (24) can be transformed to 

A^(^o) - ^ (I)' (l + 10-278 P^/^a.) ■ (25) 

Atomic interactions increase the condensate fluctuations in the trapped gas, leaving the fluctua- 
tions normal. 

Considering fluctuations in interacting systems, in the frame of the Bogolubov theory [70-72] , it 
is necessary, as is emphasized above, to omit in the course of calculations, all terms containing the 
products of more than two operators ijj, or a^, corresponding to noncondensed particles. Retaining 
such higher-order terms would remove the consideration outside the region of applicability of the 
Bogolubov approach. For example, the dispersion (22), derived in the frame of this approximation, 
is normal. However, when calculating the dispersion A'^{Nout), one deals with the fourth-order 
product of the operators a^, corresponding to the particles excited out of the condensate. The 
direct calculation of such products, with the standard replacement of the sums over momenta by 
integration in momentum space, yields infrared divergencies for both uniform as well as trapped 
gases. These actual divergencies can be formally avoided by considering the discretized sums over 
the low-energy phonon modes [76], which, anyway, results in anomalous fluctuations A'^{Nout) ^ 
N^/^ for all interacting systems, trapped or untrapped. This anomalous behaviour arises from 
the occurrence of infrared divergencies due to the term containing the fourth-order product of 
operators of noncondensed particles [76-78]. 

It is, by the way, easy to get the N'^^'^ behaviour of /S^^Nout) without discretizing the low- 
energy phonon modes. Really, the infrared divergence arising from the improper usage of the 
fourth-order operator term, which should be omitted in the Bogolubov theory, is caused by the 
integral N J dp/p^. For the lower limit of the latter, we could take Pmin = Ti/R, with the system 
radius R = {dV/AnY^^. Ten we would immediately get A'^{Nout) ~ N"^^^. 

The operator of the total number of particles is the sum N = Nq + N^ut- Therefore 

A2(7V) = A2(7Vo) + A^Nout) + 2 < NoNout > -^N^N^t ■ 

Treating the condensate field as a nonoperator quantity, one has Nq — Nq =< Nq >. Then, in the 
grand canonical ensemble, A^(7V) = A'^{Nout)- Hence, the anomalous behaviour of the dispersion 
A^(Nout) implies the same anomalous behaviour of A^(7V). The latter would mean that, in the 
thermodynamic limit, the sound velocity becomes zero, compressibility diverges, kt oo, and the 
structural factor S{0) — > oo also becomes infinite everywhere below the transition temperature. 
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Such anomalies have never been observed neither for hquid hehum nor for trapped gases. Thus, the 
excitation spectrum of the trapped Bose-Einstein condensate of ®^Rb atoms has been thoroughly 
investigated, including its long- wavelength limit [79]. The excitation spectrum, sound velocity, and 
the structural factor were found to be in perfect agreement with the Bogolubov theory predictions. 
And no anomalies of these quantities were observed. In this way, we have to conclude that the 
anomalous fluctuations in interacting Bose systems, with A'^{Nout) ~ A^^/^, are just artifacts 
caused by a not self-consistent calculational procedure. 

The problem of the fictitious appearance of anomalous fluctuations is rather general and oc- 
curs for any system with continuous symmetry, when calculations are not self-consistent. Let us 
consider a system with a Hamiltonian containing the vector operators S, whose rotation leaves the 
Hamiltonian unchanged. For concreteness, we may keep in mind the spin operators. The order 
parameter is M =< S >. An external field h acts on the operator S. The influence of the exter- 
nal field on the order parameter is characterized by the tensor of susceptibility Xap = ^Ma^/dhfj, 
whose general form can be presented as 

Xa/3 = riaUpXll + (Sap " nanp)x± , 

where n = h/h and h = |h|. Here X|| is the longitudinal susceptibility, while x± is the transverse 
susceptibility. Above the critical temperature Tc, the system is in symmetric phase, so that there 
are no distinct directions of the order parameter, which implies that x\\ = X± for T > Tc. 

But below T(., the longitudinal and transverse susceptibilities are, in general, different [80]. 
Usually, the transverse susceptibility, in a model with continuous symmetry, is larger than the 
longitudinal susceptibility, which results in the dominance of directional over size fluctuations in 
destroying order [81]. It may even happen that the transverse susceptibility diverges as x± ~ h~^, 
when h ^ 0. This occurs because, under continuous symmetry, all directions of the vector 
order parameter are equivalent, and an infinitesimally small transverse field may turn the order 
parameter by a finite angle. Infinite susceptibility means that the considered model is unstable. 
Such an instability can be easily removed by adding to the model a finite additional field breaking 
the continuous symmetry. Actually, any real system always has such breaking-symmetry field. For 
example, the isotropic Heisenberg model is well known to be a cartoon of real magnetic systems, 
always possessing anisotropic terms due to spin-spin dipolc interactions, spin-orbital interactions, 
indirect exchange, crystalline fields of different nature, and demagnetization anisotropic factors 
caused by the shape of the sample [82-84] . In the presence of a finite, though may be weak, field 
h, the transverse susceptibility x± ~ is finite. The related transverse correlation function 
C_L =< SS±SS± >, where SS± = S±— < S± >, has the spatial dependence 

'^^(^) ^ (r ^ OO) , 

with the correlation length ^ ~ h~^. The corresponding Fourier transform reads C±{q) ~ (q^ + 
as 5 — > 0. At zero field h = 0, the form C±{q) ~ 1/g^ is known as the Goldstone theorem, 
whose rigorous proof for Bose and Fermi systems was given by Bogolubov [72]. 

For the longitudinal correlation function C\\ =< ^^ii^^*!! >, with = 'S'||— < S\\ >, one may 
write [80] an approximate form C\\ ~< {SS±y{SS±y >. Treating the latter in the hydrodynamic 
approximation, one gets [80] the Wick- type decouphng 
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This gives the relation C\\{r) ~ Cj_{r), from which it follows that, at zero field h = 0, one 
has C\\{r) ~ l / r'^i'^-'^) and C||(g) ~ l/q^^'^. Recall that the index 77 = in the hydrodynamic 
approximation. This tells that the longitudinal susceptibility x\\ ~ h'^^"^ also diverges at zero 
field [80]. The character of this divergence follows from the relation x\\ ~ / C\\{r)dr ~ N^'^-'^)/^^ 
valid for 2 < d < 4. Then the dispersion of the spin operator S, for which A^(S) ~ -^X||, behaves 
as A2(S) ~ iV('^+^)/3. In the three-dimensional space, one gets A'^(S) N^/^ for all T < T^, that 
is, the same kind of anomalous fluctuations as for Bose-Einstein condensates. 

If such anomalous fluctuations would really exist, this would mean that there arc no stable 
ordered phases for the systems with broken continuous symmetry. This kind of conclusion would 
certainly be incorrect. It is well known that stable ordered phases, such as magnetic phases, do 
exist [82-85]. The longitudinal susceptibility for these i^liascs diverges solely at the critical point 
Tc, where the correlation length behaves as ^ ~ a|l — T/T^^'^, with a critical exponent u, and a 
being the mean intcrparticle distance. Far below T^, the correlation length at zero field is always 
finite, the correlation function exponentially decays, and the longitudinal susceptibility is also 
finite, hence, the order-parameter fluctuations are normal. 

The origin of the fictitious anomalous fluctuations in the systems with continuous symmetry 
is the same as that occurring in the case of the Bose-Einstein condensate. The hydrodynamic 
approximation consists in leaving in the Hamiltonian only the terms up to second order in the 
operators 6S±. For spin systems, ^^j, plays the role of the magnon operators. To calculate A^(S), 
one needs to consider the product of four operators SS±, which one treats by invoking the Wick de- 
coupling and using the hydrodynamic approximation. However, the four-order operator products 
cannot be correctly treated in the second-order hydrodynamic approximation. This is the same as 
to try to employ a linear approximation for describing strongly nonlinear effects. Because of this 
inconsistence, there arise infrared divergencies leading to the appearance of fictitious anomalous 
fluctuations. 

The same inconsistence appears if one tries to calculate the four-order terms in the parti- 
cle dispersion of quantum liquids, such as helium, by using the hydrodynamic approach, which 
is a second-order approximation. No anomalous fluctuations arise, when calculations are self- 
consistent, as is demonstrated by Bogolubov [71] when considering the fluctuations in weakly 
interacting gases with Bose-Einstein condensate. 

Let us now summarize the behaviour of fluctuations in Bose-condensed systems. The ideal 
uniform Bosc gas, below the condensation temperature, is an unstable system with anomalous 
condensate fluctuations described by the dispersion (20). Switching on atomic interactions stabi- 
lizes the system making fluctuations normal according to Eq. (22). 

The trapped ideal Bose gas is stable, with normal fluctuations, as is described by the dispersion 
(24) for a harmonic confining potential. Switching on interactions leaves the fluctuations normal, 
but increases them, as is illustrated by Eq. (25). Different properties of ideal uniform and trapped 
gases explain why the condensation temperature in a uniform gas displays the reentrant behaviour 
as s function of the scattering length a^, while the condensation temperature in a trapped gas is 
a monotonely decreasing function of a^. 
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4 Condensate Growth 



Formation of the Bose-Einstein condensed phase from a strongly nonequihbrium, completely dis- 
ordered, initial state, ultimately leading to long-range ordering, is one of the most interesting and 
fundamental dynamic problems of the physics of multiparticle systems. The attempts of describing 
this process have encountered controversies. First, Levich and Yakhot [86], using the semiclassical 
kinetic equations, concluded that the process of condensation requires infinite time. To make this 
time finite, it is necessary to assume the existence on the nuclei of the condensed phase in the be- 
ginning of the cooling process. Svistunov [87] and Kagan et al. [88-90] also followed a semiclassical 
approach based on the nonlinear Schrodinger equation. Their analysis of the kinetics starts with 
the statement that the weakly interacting Bose gas possesses large, relative to unity, occupation 
numbers and, thus, is essentially a classical field phenomenon. According to their picture [87-90] , 
the condensate formation goes through three stages. At the first stage, substantial suppression of 
nonequihbrium fluctuations of the density occurs. Then a kind of short-or medium-range order 
appears, when the quasicondensate is formed, such that the density fluctuations are suppressed, 
but phase fluctuations arc essential. And the true condensate arises only at the last, third, stage, 
when the long-range order develops, with all fluctuations of density as well as phase being relaxed. 
Stoof [91,92] stressed the importance of quantum fluctuations at the initial stage. 

Bose-Einstein condensate is a coherent system. The development of coherence in a nonequi- 
hbrium Bose gas is somewhat analogous to the emergence of coherence in radiating laser-type 
systems, where the coherent stage is preceded by the quantum stage, with quantum fluctuations 
triggering the appearance of the coherent component [93,94]. 

In order to better understand how the condensate is being formed, it is important to recollect 
the main characteristic length and time scales that are common for all many-particle systems 



One evident characteristic length is the interaction radius, which can be well represented by a 
scattering length a^. This defines the effective particle velocity v ~ h/mas, with m being particle 
mass. The interaction time 



describes the time during which two particles strongly interact with each other in the course of 
mutual scattering. The related interval 



is called the dynamic stage, when particles move more or less independently from each other, their 
motion being described by dynamical equations. 

The mean free path A = l/po^, where p ~ a"^, with a being the mean interatomic distance, is 
the characteristic length between two scattering events. Typical relations are 



[95-97]. 



<t < Tint 



{dynamic stage) 



(27) 
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The correlation time or local- equilibrium time is 



A 

Tloc = - 



V hpas 
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(28) 
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It is connected with the interaction time (26) as 



'^int 3 
Tloc 



The local-equihbrium time (28) can also be expressed through other characteristic quantities. The 
heahng or coherence length ^ = h/ms is defined through the sound velocity s ~ hy/pa^/m. The 

latter is related to the particle velocity as s/v ~ y^paf. For the coherence length, one has 

=, > T ~ V ■ 

ypdg a y Os A 

The typical energy scale is given by £ ~ fi^pas/m. Therefore, for the local-equilibrium time (28) 
we get 

\ ^ h 
V s e 
The temporal interval 

'Tint < t < Tioc {kinetic stage) (29) 

is the kinetic stage, during which particle motion can be described by kinetic equations. At this 
stage, because of mutual interactions, interparticle correlations begin arising. 

The existence of the next stage was advanced in [98] (see also [95-97]). Then the overall system 
consists of parts inside each of which there is some kind of ordering, while different parts are not 
correlated with each other. The typical size of each part, //, is mesoscopic, being in between the 
mean particle distance a and the linear size of the whole system, L, so that 



T/oc = - ~ - ~ 



a^lf^L. (30) 

Since these mesoscopic subregions are not correlated with each other, but distributed randomly 
in space, they are termed heterophase fluctuations. Examples of the latter are magnetic clusters 
inside a paramagnetic matrix, disordered regions inside a crystal, crystalline formations inside a 
liquid, superconducting droplets inside a normal metal, subregions with differing local densities, 
and so on [95]. In particular, these could be separated spatial regions, inside each of which there is 
a coherent Bose-Einstein condensate, but with no coherence between the different spatial regions. 
Denoting the average velocity of motion for each heterophase nucleus as Vf, we may introduce the 
heterophase fluctuation time 

r, = ^, (31) 
which is a characteristic lifetime of a local heterophase fluctuation. During the temporal interval 

'Tioc <^t<^Tf {heterophase stage) , (32) 

defining the heterophase stage, the system is spatially nonuniform, consisting of subregions with 
different order parameters [95] . The actual value of tj depends on the particular system considered, 
but it is always larger than the local-equilibrium time tioc, since to organize a kind of order across 
the distance If ^ a requires a time t/ ^ tjoc- 
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The final temporal interval is the hydrodynamic stage lasting the time defined by the inequal- 
ities 

Tf <^t <^ Texp {hydrodynamic stage) , (33) 

where Texp is the experiment time or the observation time. In many cases, one may set r^xp oo. 
At this final stage, nonequilibrium systems can be described by hydrodynamic equations. 

For usual condense matter, such as solid or liquid, one has pa^ ~ 1, so that s v and 
Os ~ a ~ A ~ ^ <^ Z/. Then the kinetic stage (29) is practically absent, and only three stages are 
left, dynamic, heterophase, and hydrodynamic. 

For dilute gases, when pa^ <^ 1, one has s -C so that <^ a <^ ^ -C A. Then all 
four temporal stages are present, since Tint ^ tioc t/. To estimate the characteristic scales 
for dilute gases, let us take the values typical of experiments with the trapped atoms of ^^Rb 
and ^^Na, for which m ~ IQ-^^ ^ 5 x 10"^ cm, and p ~ (10^^ - 10^^) sm-^. Then the 

parameter pa^ ~ 10~^ — 10~^ is really small. Atomic velocity is f ~ 10 cm/s and sound velocity 
is s ~ (10~^ — 10^^) cm/s. Mean interatomic distance is a, ~ (10^^ — 10^"^) cm. Coherence length 
is ^ ~ (10~^ — 10^^) cm. Mean free path is A ~ (10^^ — 10^) cm. The interaction time (26) is 
Tint ~ 10^^ s, that is, the dynamic stage (27) is very short. The local- equilibrium time (28) is 
Tioc ~ (10~^ — 1) s. Hence, the kinetic stage (29) is much longer than the dynamic stage (27), 
because of which the latter can be neglected in practical description. 

In the frame of the defined time scales, one may consider the formation of the Bose-Einstein 
condensate from a nonequilibrium disordered state. Then, according to [87-90], we have the 
following picture. The dynamic stage (27), being rather short, is usually omitted from the con- 
sideration. During the kinetic stage (29), substantial suppression of nonequilibrium fluctuations 
of the density occurs. After the local-equilibrium time tioc, mid-range order develops in regions of 
size If > ^. Different spatial regions are not correlated with each other. Such a situation is typical 
of heterophase fluctuations. At the heterophase stage (32), density fluctuations are suppressed, 
but phase fluctuations are essential. The authors [87-90] call such a system quasicondensate. Af- 
ter the fluctuation time t/, coherence penetrates through the whole system, true long-range order 
develops, so that the Bose-Einstein condensate is flnally formed. At the last hydrodynamic stage 
(33), the Bose-condensed gas relaxes to equilibrium. 

The formation process of a Bose-Einstein condensate in a trap was described by Gardiner et 
al. [99,100] using a master equation based on quantum kinetic theory. The description assumed 
the existence of one condensate mode in interaction with an equilibrium thermalized bath of 
noncondensate atoms at flxed temperature T and chemical potential p. The number of condensed 
atoms A^o(^) satisfles the rate equation 

No] , (34) 

in which Ap = Pc{Nq) — p is the difference between the condensate chemical potential Pc{No) 
and the bath chemical potential p, while W{No) is a transition probabihty from the bath to the 
condensate. A collision between a pair of atoms initially in the bath of atomic vapor results in one 
of the atoms going to the condensate, because of which W{0) is not zero. Thus, with the initial 
condition iVo(O) = 0, at the initial stage, Eq. (34) gives iVo(t) ~ 2W{0)t. 

Growth of a Bose-Einstein condensate from thermal vapor was experimentally studied [101] 
for atoms of ^^Rb. It was found that there exists a latency time of (10"^ — 10~^) s at which 



dt 



2W{No) \ 1 + 
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a condensate is first detected. According to the time scales described above, the latency time 
should correspond to the local-equilibrium time r^oc, since before this time, at the kinetic stage, 
there cannot be yet any condensate, which appears after r|„,.. Then the slow linear growth of the 
condensate was observed [101] till the time of 5 x 10~^ s. The latter should correspond to the 
heterophase fluctuation time t/. The interval Tioc < t < Tf is the heterophase quasicondensate 
stage. The numerical calculations based on the rate equation (34) do not reproduce the measured 
condensate growth data accurately in the quasicondensate regime. This equation describes better 
the condensate growth at the last hydrodynamic stage, when t ^ tj. Then there is exponential 
growth till the condensate reaches its asymptotic equilibrium value at about 10^^ s. At this last 
stage, instead of Eq. (34), one can use the simple form 

N^{t) = 7Vo(io) {1 - exp [-V{t - t^)]} , (35) 

valid for ^ t/oc, with a fitted relaxation rate F. 

There is no yet a quantitative theory describing noncquilibrium condensate formation from 
the very beginning of the process, starting from the kinetic stage, with no condensate, and going 
through the heterophase and hydrodynamic stages. 



5 Low-Dimensional Condensates 

The properties of systems in one and two dimensions can essentially differ from those in three 
dimensions. The Bose systems with the contact interaction 

$rf(r) = Aa5{r) , (36) 

where r is a (i-dimensional vector, have been studied for the spatial dimensionality d = 1 and 
(i = 2 as well as for the usual three-dimensional case. For the latter, one has As = Anh^as/m. 

First of all let us recall that the Bose- Einstein condensation in uniform systems does not exist 
at finite temperatures in d = 1 as well as in o? = 2. This rigorously follows from the Bogolubov 
theorem [72] as is shown in [102,103]. 

The features of low- dimensional Bose systems with the interaction (36) can differ owing to 
the strength of their atomic interactions. This can be illustrated by the classical example of the 
one-dimensional Bose gas, solved by Lieb and Liniger [104,105]. Let us consider the ground-state 
energy E, written in the dimensionless form as 

being a function of the dimensionless coupling 

mAi 
n Hi 

In the weak-coupling limit, when g ^ I, the energy e{g) possesses an asymptotic expansion 

e{g) c^g + Cg/ac/^/^ + C2^^ + C5/2g^/^ , 



20 



with the coefficients [106-108] 



4 

C3/2 = = -0.424413 , C2 = 0.06535 , C5/2 = -0.017201 . 

3n 

In the strong couphng hmit, with g —>■ 00, which is also called the Tonks-Girardeau regime [109- 
111], one has 

e(oo) = — = 3.289868 . 

^ ' 3 

In the Tonks-Girardeau regime, the wave functions of the Bose system are exactly mapped [110] 
by the Girardeau mapping = where are the eigenstates of the ideal one-dimensional 
Fermi gas provided by the simple Slater determinants. Thermodynamic functions of the one- 
dimensional system are analytic at any finite temperature [112], in agreement with the absence of 
a finite-temperature phase transition [102,103]. 

Though there is no Bose-Einstcin condensates in d = 1,2 in finitc-tcmpcraturc uniform sys- 
tems, the properties of the latter display noticeable change when varying temperature. This can 
be demonstrated by considering the first-order density matrix 

pS'^)(r)^<V^t(r)V;(0)>, (37) 

in which the upper index d shows dimensionality. 

For the one-dimensional case, at zero temperature, one finds [113-115] the large-distance decay 

pi^)(r)~r-^ (T = 0), (38) 

where 7 = mc/27rp and c = ^J^A^Jm is the sound velocity. While at any finite temperature 
[115,116], for r — > 00, one has 

p(^)(r)~exp(-!^r) (T > 0) . (39) 

This means that there is no long-range order at any T > 0, but at T = there arises mid-range 
order characterized by the power-law decay (38). Therefore, at T = 0, there can exist superfiuidity. 

For the two-dimensional uniform gas at zero temperature the long-range order does exist [5], 
since 

pf)(r)^Po>0 (r = 0), (40) 

as r — i> 00. At finite temperatures, below the Kosterlitz-Thouless temperature Tkt, which is of 
the order of the degeneracy temperature ~ h^p'^^'^/mkB, the matrix (37) decays algebraically 
[117-120], 

pf)(r)~r-^ {0<T<Tkt), (41) 

similarly to the case (38), where r — > 00. Above the temperature Tkt, the density matrix (37) 
decays exponentially, in analogy with Eq. (39). Therefore, in the two-dimensional uniform Bose 
systems, there is Bose-Einstein condensate at T = 0, and there exists superfiuidity at < T < Tkt, 
though there is no condensate at any finite temperature. 

The situation is difi'erent for low-dimensional trapped gases, for which Bose condensation can 
occur at finite temperatures. Recall that for a finite number of particles there is never a rigorous 
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phase transition, but Bose condensation happens as a gradual crossover, which can be very sharp 
and practically undistinguishable from the rigorous phase transition, when the number of particles 
N is sufficiently large. This concerns trapped systems of any finite dimensionality. 

Phase transformations of the crossover type can occur, in some cases, even in the thermo- 
dynamic limit [95,121]. Therefore, it is useful to generalize the notion of the phase transition, 
according to the behaviour of the order parameter, as follows [95]. The phase transition is of first 
order, when the order parameter changes, at the transition point, by a jump. The phase transition 
is of second order, when the order parameter at the transition point changes continuously from 
zero to a finite value. And the phase transition is of third order, or of crossover type, when the 
order parameter changes continuously, being nonzero, so that there occurs a qualitative change of 
its behaviour. The crossover point can be defined as a point separating the qualitatively different 
regions. Often, the crossover points can be identified with the locations of the extrema of some 
derivatives [95,121,122]. 

Bose-Einstein condensation of trapped ideal gases in different dimensions and in various traps 
has been reviewed in [2]. For d — 1, the quasiclassical approximation does not show any condensa- 
tion [2,123]. However, treating accurately the lower discrete levels displays a crossover transition, 
which, for a harmonic trap of frequency cuz, can be located [124] at the temperature 

(1) _ hu^N 

For a two-dimensional harmonic trap, with the radial frequency uj±, the quasiclassical approach 
is appropriate yielding [2,123] the crossover temperature 



TTfcs 

Below the temperatures (42) or (43) there occurs a macroscopic occupation of the single-particle 
ground-state energy level, because of which it is permissible to talk on Bose-Einstein condensation 
[2]. The crossover temperature in low-dimensional trapped dilute gases with interaction is of 
the order of the temperatures (42) and (43). 

Dealing with interacting gases, it is not always easy to define the single-particle spectrum, 
from which one would like to separate out the groimd-state level ascribed to Bose condensate [2]. 
Nevertheless, in many cases, the condensate density can be defined, serving as a reliable order 
parameter. A more general approach could be by introducing the order indices [65,125]. An order 
index, for a reduced p-particle density matrix pp, is defined [125] as 

where || ■ || implies the Hermitian norm. If the total number of particles N is large, one has 
In |Trpp| ~ pin A^. For a pure Bose condensate, Up = 1, while uUp = when there is no condensate. 
And in the intermediate region, the order index (44) varies in the interval < a;^ < 1, being equal 
to LOp = In No/ In N [125]. 

Although in reality there are no exactly one-or two-dimensional systems, it is possible, be 
varying the shape of the trapping potentials, to produce dilute cold gases in highly anisotropic 
configurations, where the motion of atoms is quenched in one or two directions. The radial 
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frequency a;_L and the axial frequency ujz of a cylindrical harmonic trap define, respectively, the 
transverse and longitudinal oscillator lengths 



U 



! h 

muj± 



L 



When the trap is strongly squeezed in the radial direction, so that l± <^ Iz, one has the quasi- 
one-dimensional configuration. And if the trap is squeezed in the axial direction, so that Iz <^ l±, 
one gets the quasi-two-dimensional system. These two opposite configurations correspond to the 
cigar-shape or disc-shape traps. For such highly anisotropic traps, it is possible to reduce the 
generic three-dimensional consideration to an effective one- or two-dimensional case, with the 
appropriate interaction constant in the potential (36). If the scattering length Ug is smaller, by 
modulus, than all oscillator lengths, the reduction to a low-dimensional case is straightforward, 
requiring just to integrate out the variables related to the quenched directions [2]. Thus, for the 
quasi-one-dimensional trap, one finds 



A, 



(|a,| < < Iz) . 



And for the quasi-two-dimensional trap, one obtains 



mL 



(|a,| < < l±) . 



(45) 



(46) 



The situation becomes more complicated when the scattering length is comparable with an 
oscillator length. For example, if the scattering length is of the order or even larger than l±, then in 
the quasi-onc-dimcnsional case, the quartic self- interaction could be replaced by a nonpolynomial 
interaction [126-128] or by a sextic term [129,130]. A simple form of Ai for arbitrary scattering 
lengths was suggested by Olshanii [131], 



2h^as 



IC(l/2)|a. 
V2l± 



(47) 



where C(l/2) — —1.46. In the limit \as\ <^ l±, this reduces to Eq. (45). Bose gases confined in 
highly elongated harmonic traps were investigated over a wide range of interaction strengths using 
quantum Monte Carlo techniques [132]. The properties of the Bose gases under tight transverse 
confinement were found to be well reproduced by a one-dimensional Hamiltonian with the contact 
interaction (36) and the interaction parameter (47). The latter well represents atomic interactions 
for arbitrary scattering lengths, positive as well as negative. For |as| » lj_, the parameter (47) 
tends to the unitary limit 



A, 
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which does not depend on and is attractive. 

In the quasi-two-dimensional case, the effective interaction parameter for arbitrary ag was 
considered by Petrov et al. [133,134]. Their results can be presented in the form 
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(49) 
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In the limit of the asymptotically small scattering length \as\ <^ Iz, Eq. (49) reduces to the form 
(46). However, under the condition 



the interaction parameter (49) becomes 

A. ^ , (50) 

mln(/_L/ vStt/z) 

and does not depend on a^,. One could also introduce, instead of the contact interaction (36), 
some more complicated interaction potentials [130,135,136]. 

Note that for repulsive interactions, when > 0, the parameter (49) is always positive. But 
for attractive interactions, when < 0, the effective quasi-two-dimensional interaction (49) has 
a resonance form, diverging at an intermediate scattering length 

a. — 



2\n{ljV27rlz) ' 

and changing sigh from negative to positive. 

Quasi-one-dimensional strongly alongated cigar-shaped traps have been used in several experi- 
ments [137-145]. Quasi- two-dimensional disc-shaped traps have also been employed in experiments 

[141,146]. For example, in experiments [142,143] the trap frequencies were uj± = 271 x 365 Hz and 
uJz = 271 X 14 Hz. The frequency ratio uo^/uOz = 26 corresponds to the length ratio lz/l'± = 5, which 
describes an alongated trap. The radial squeezing can be increased [139], reaching u!± = 27t x 715 
Hz, which translates into the length ratio lz/l± — 7. Quasi-one-dimensional condensates can also 
exist in waveguides [147] . 



6 Crossover Region 

In the vicinity of critical points, there always exists a region, where the system experiences strong 
fluctuations. This region is termed critical. In finite systems, the phase transition is of crossover 
type. Around the point of crossover, there can arise strong fluctuations, similarly to those ap- 
pearing in the vicinity of the critical points. Then the region around a crossover point, where 
fluctuations are strong, can be called crossover region. Bose-Einstein condensation of trapped 
atoms is always a crossover. In three dimensions, the latter can be rather sharp, when the num- 
ber of trapped particles is large. Therefore the crossover region can be quite narrow. But in 
very anisotropic traps, in which the quasi-one-dimensional or quasi-two-dimensional configuration 
is realized, the condensation crossover can go as a very smeared process, with a wide crossover 
region. 

Considering phase transitions under varying temperature, one usually defines the critical region 
by means of the Ginzburg temperature Tg, such that in the critical region Tq <T <Tc fluctuations 
are strong, while essentially below Tg, they are suppressed [148]. In general, one may distinguish 
the longitudinal and transverse fluctuations of an order parameter. The transverse, or directional, 
fluctuations often survive much longer, when lowering temperature below Tc, then the longitudinal, 
or size, fluctuations. Therefore the actual critical region is characterized rather by the transverse 
then longitudinal fluctuations [81]. 
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In Bose- Einstein condensates, one may also distinguish two types of fluctuations. The size fluc- 
tuations of the order parameter correspond to density fluctuations. The transverse fluctuations 
are usually associated with phase fluctuations. In the vicinity of the condensation temperature 
Tc, both these fluctuations are strong. Lowering temperature below T^. flrst diminishes the density 
fluctuations, while the phase fluctuations remain yet strong. Lowering temperature further, one 
comes to a point Tj, where phase fluctuations are getting suppressed as well. Thus, the crossover 
region is deflned by the temperature interval Tf < T < T^. For temperatures essentially lower 
than Tc, but yet much higher than T/-, the system can be described as being spatially divided into 
subvolumes, each of which being a coherent condensate, but with no coherence between different 
subvolumes. Such a system is an example of heterophase matter [95,149], in this particular case, 
of a heterophase condensate. It has also been called quasicondensate [150]. As is discussed in Sec- 
tion 4, the heterophase condensate, or quasicondensate, appears in the process of nonequilibrium 
condensation. And it may also exist in the crossover region of equilibrium systems. The crossover 
regions are especially pronounced in low-dimensional trapped gases [133,151,152], though they 
also occur in purely three-dimensional cases [153]. 

To consider phase fluctuations, it is necessary, first of all, to introduce a phase operator. In 
the works on Bose-Einstein condensation, this is commonly done by presenting the field operator 
■0 as 

^/'(r) = exp{i^(r)} yp(^ ? (51) 

Here (p is the phase operator and p = ■ip'^il). The question-mark in Eq. (51) is the sign of warning 
that one should not blindly trust this equation, which, strictly speaking, is not correct. The 
problem of introducing the phase operator is a long-standing and not yet completely solved one, 
as can be inferred from reviews [154-157] (see also a recent discussion in [158]). Thus, in order to 
satisfy the identity p = the phase operator is to be Hermitian, (p^ — (p^ so that the creation 
field operator could be written as 

^\v)^{^)^^v{-i(p{v)) ? (52) 

Here and in what follows, we continue marking by the question-mark those equations that in strict 
sense are not correct. Above all, the phase operator cannot be defined as a Hermitian operator 
[154-157]. 

The fact that the representation (51) is not correct can be demonstrated as follows. From this 
representation, one has exp(i<^) = ■\\)p~^l'^ . Then, using the commutation relation 

[V^(r), p(r')] = V^(r)5(r - r') , 

which is vahd for any (either Bose or Fermi) field operators, one gets 

[e^^W, p(r')] = e'^W5(r - r') ? (53) 

From here, after expanding the exponentials, it follows the commutation relation 

[p(r), (p{v')\=%b{j-v') ? (54) 

This is what one standardly employs announcing the density and phase operators being canonically 
conjugated. However, the relation (54) is improper. Really, if it were correct, then for the number- 
of-particle operator N = J p{r)dr one would have 

[N, ip{r)] = I ? (55) 
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The latter equation is already suspicious, with its left-hand side depending on r and the right- 
hand side containing no such a dependence. Moreover, Eq. (55) can be reduced to a completely 
ugly form by taking its matrix element over the states \n > from the number basis, for which 
N\n >= n\n >. This yields 

(n — n') < n\(p\n' > — i5nn' ? (56) 

Setting here n — n' results in a senseless equation — i. 

In this way, the representation (51) cannot possess a rigorous mathematical meaning. Then 
why people use it? To answer the question, let us note that if, instead of the field operator ■0, we 
would deal with a classical wave function, then the form ip = y^e*''' would be absolutely proper, 
but with the density p and phase (p being also classical functions. This remark suggests that 
the representation (51) could be accepted as an approximate form for a kind of semiclassical ap- 
proximation, when the density and phase operators could be treated as almost classical functions. 
This is equivalent to saying that the density and phase fluctuations above their average values, 
for which Sp = p— < p > and S^p = 0— < if >, should be very small. The situation characterized 
by small fluctuations corresponds to the hydrodynamic approximation. 

There could be several equivalent ways of introducing a semiclassical picture. One could treat 
the evolution equations for the mode operators as semiclassical, in analogy with the semiclassical 
laser theory [159]. One could start from the Gross-Pitaevskii equation, present the wave function 
in the form ip = y/pe*''', and then quantize the classical density and phase [4,160]. Or one could 
employ from the very beginning the representation (51), but not forgetting that it is only ap- 
proximately meaningful in the hydrodynamic limit, when all fluctuations are small. In this way, 
one often treats the density and phase fluctuations as not correlated with each other, permitting 
their separate averaging. And the phase operators are usually considered as random Gaussian 
variables, which allows one to invoke, for a function f{if) of ip, the equality 



< exp{f{0)} > = exp |i < f{ip) >| . 



It is worthwhile to emphasize that the phase, as such, is not an observable quantity. What one 
is able to observe, except the density of atoms, is the density of current, whose operator reads 



I _ ih 
2m 



or the velocity, for which the velocity operator is v = j/p. For the classical quantities, one would 
have j = {hp/m)'Vip and v = {h/m)'Vip. Respectively, instead of talking about unobscrvable 
phase fluctuations, one could discuss the current or velocity fluctuations, which can be correctly 
deflned without mentioning the problematic phase operator. 

The density operator p and the operator of current j are connected by the continuity equation 

dp , ^ 

which is valid for an arbitrary trapping potential and any atomic interactions. Hence, current 
fluctuations (or phase fluctuations) would, generally, provoke density fluctuations, and vice versa, 
since they are also connected by the continuity equation 

d 

—6p + div5j = . 
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Density fluctuations always excite current fluctuations. However, there can exist current fluctua- 
tions that do not initiate density fluctuations. This takes place for the divergenceless fluctuations, 
for which div^j = 0. Hence, there can exist circulating current fluctuations that leave the atomic 
density undisturbed. This explains why the crossover region is mainly defined by the phase (cur- 
rent) fluctuations, but not by the density fluctuations. 

In terms of the observable quantities, we can say that, in the vicinity of the crossover temper- 
ature Tc, both the density and current fluctuations are strong. In the part of the crossover region, 
where Tf < T <^ Tc, the density fluctuations can be suppressed, but the current fluctuations are 
yet essential, which is typical of quasicondensate. And for T <^ T/, true condensate develops. 

For quasi-one-dimensional gases, the fluctuations temperature 7/ can be estimated [151,152] 

as 

^ Kb \aslzj 
The characteristic phase-fluctuation length is 

/;2 ;2\ 1/3 

(58) 




which should be much less than the length L of the trap. 

In quasi-two-dimensional gases, the fluctuation temperature becomes [133] 




(59) 

ATexp(A|p'/') , (60) 



with the fluctuation length 



where it is assumed that <^ 1^ and At = \j2nh /mkBT is the thermal wavelength. Here, it 

should be that ly <^ R, with R being the trap radius. 

The existence of fluctuations in the crossover region of elongated Bose-Einstein condensates 
has been confirmed experimentally [142-145]. 



7 Condensates and Superfluids 

The relation between Bose-Einstein condensation and superfiuidity is a long-standing and not 
completely understood problem. One often thinks of these two phenomena as being closely related. 
A detailed description of the historical evolution of views on this problem has recently been done 
in [2]. A brief and nice discussion can also be found in [161]. The current understanding [2,161] 
is that Bose-Einstein condensation is neither necessary nor sufficient for superfiuidity. The latter 
can occur without Bose-Einstein condensation. Thus, in two-dimensional uniform systems, as is 
mentioned in Section 5, there can exist superfluidity, though the condensate cannot appear at 
any flnite temperature. Huang [161] gives an example of a Bose system in a random fleld, when 
superfluidity is destroyed, even though there is a Bose condensate. Carusotto and Castin [162] 
find superfiuid behaviour in a one-dimensional ring of weakly interacting gas outside the Bose 
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condensed regime. In the present section, we shall adduce arguments why there can be no simple 
and general relation between the densities of condensate and of superfluid. 
The condensate density 

PO = P - Pcmt (61) 

is defined as the difference between the total density and that of atoms outside the condensate. 
The latter density for a uniform gas reads 

P<mt = TT^^k , rife = < alttk > . (62) 

For instance, in the Bogolubov approximation at low temperature one finds [70,71,163] 

Po 8 /-^ miksTf 

— — 1 - -;r-F= \ po-i - ^ , (63) 



where c = y Anh^pas/m is the Bogolubov sound velocity. 

To find the superfluid density Ps, one consider a system boosted with the velocity v = {v"}. 
In the frame, moving with the same velocity, the Hamiltonian is a functional H{il)} of the field 
operators ■0. In the laboratory frame, the Hamiltonian Hy — H{ipy} is the same functional, but 
of the field operators ipy given by the Galilean transformation 

V'^(r, t) = ip{r, t) exp | ^ (^mv ■ r - ^ | . 

The superfluid component is interpreted as that part of the system, which nontrivially responds 
to the velocity boost, with the superfluid density 

,.^;^]iS,lEi,<K>,. (64) 
where the averaging, with the Hamiltonian Hy, is over the total momentum 

Pv = J ijl{r){-thV)ijy{r) dr. 
Performing the limit i> — > in Eq. (64), one has 

Ps = P- Pn , (65) 

with the normal density 

- £n ^"f*") (^ - i^) 

expressed through the momentum dispersion 

A2(P) =< P2 > - < P >2 . 

When the total current across the system is absent, then < P >= 0. Combining Eqs. (61) and 
(65), one gets the formal relation 

Ps = P0 + Pout - Pn ■ (67) 
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Note that the derivative in Eq. (64) is a particular case of differentiation with respect to a 
parameter. In general, if a parameter A enters an operator A — A{X) as well as the Hamiltonian 
H = H{X), then 

d<A> dA ^( J dH \ 



d\ dX ^ \ dX dX 



This formula of parametric differentiation is useful for different particular applications. 

To find the normal density (66) one needs to calculate the momentum dispersion, which, for a 
motionless gas, is A^(P) =< >. This can be expressed [2] as 

<P'^ > = h'^ lim lim V3 • V4 [pi(r4, r3)5(ri - rg) - P2(r3, r4, ri, rg)] dridr2 , (68) 

where the first-and second-order density matrices, respectively, are 

Pi(r4,r3) = < '0^(r3)V'(r4) > , 

P2(r3,r4,ri,r2) = < ip^{r2)ilj^{ri)ilj{r3)ilj{r4) > . 

Thus, the definition of the superfiuid density (65) involves the second derivatives of the first-and 
second-order density matrices, while the condensate density (61) requires solely the first-order 
density matrix. For interacting systems, there is no a simple general relation between the reduced 
density matrices of different orders (sec [65] and recent discussions in [164-168]). This means that 
it is hardly probable that there could exist a general relation between the condensate density and 
superfiuid density. 

For Bose systems with condensate, one usually employs the Bogolubov shift ilj{r) = + tp{r), 
separating the condensate part ■00 from the field operator ■0(r) of noncondensed particles. The 

condensate term ipQ is treated as a classical function, which becomes a constant for uniform gases. 
The field operator of noncondensed particles, for which < ip >= 0, can be expanded over a basis 
(^fc(r), so that '^(r) = I]fc^o <^A;<^fc(i')- For uniform systems, which are considered in what follows, 
the basis is formed by plane waves 93fe(r) = e^^'^/\/V. In this way, the single-particle density 
matrix reads 

Pi(ri,r2)=po + :^E^fee^''"'"' (69) 

where ri2 = ri — r2 and the property < aj^flg >= SkqUk is taken into account. 

Constructing the second-order density matrix in the Bogolubov approach, one has to omit all 
terms containing the powers higher than those of second order in the noncondensed field operators 
ijj or Ofc with k ^ 0. Then, keeping in mind the momentum conservation, resulting in the condition 
< akQq >— < akQ^k >) one obtains 



P2(r3,r4,ri,r2) = pl+ 



+77 E {< ^ka-k > (e^^ "'' + e^^ "'^) + rik (e-^ ''^^ + e-^"*""" + e-'^ ''^ + e''^-''^) } . (70) 

The diagonal element of the density matrix (70) gives the pair correlation function 

girn) = 4 P2(ri, r2, ri, ra) . (71) 
P 
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For the latter, when No Ki N, one finds 

g{r) = 1 + - y K+ < aka-k >] e^'^-' j—^ . (72) 
Prom here, one has the equahty 



/ 



2 

[^(r) - 1] dr = - lini [nk+ < a^a^k >] , 



which, in the Bogolubov approximation [70,71], yields 




with c being the Bogolubov sound velocity. It is this equation that was obtained by Bogolubov 
when calculating the particle dispersion (22), with the use of Eq. (18). A similar way can be 
followed in the calculation of the momentum dispersion (68), which is necessary for finding the 
normal density (66). 

Note that for the ideal gas it is easy to show that po = Ps (see details in [2]). This result, 
however, is not of much meaning. First, the ideal Bose gas cannot be superfluid, since it does not 
satisfy the Landau criterion of superfluidity. Second, the ideal Bose-condensed gas, as is explained 
in Sec. 3, is an unstable system, having infinite compressibility. 

It is also easy to show that, for any system at asymptotically high temperatures, one has 
Pn — p, hence Po = Ps = 0. This immediately follows from the Maxwell-Boltzmann distribution 
law rife ~ exp{P{p — cuk)}, where, assuming the parabolic spectrum cuk, the chemical potential is 

3 mkeT 
/X ~ - - /csTln ^^-^ 



2 " 27rrp2/3 

To illustrate the way of calculating the momentum dispersion (68) for interacting gases, let us 
resort to the Bogolubov approach [70-72]. Then, the annihilation and creation operators of non- 
condensed particles in the momentum representation are and a\, with k 7^ 0. After omitting in 
the Hamiltonian all terms of third and fourth order with respect to these operators, the remaining 
quadratic form is diagonalized by means of the Bogolubov canonical transformation 

Ofc = Ukh + Vfc&lfc , 4 = Ukhl + Vkh-k , 
where the operators correspond to excitations. As a result of this diagonalization, one gets 

1 

Uk — I , Vk — I , 

where 



Sk-^el + m^c'^ 
Jjk — 



and £fe is the Bogolubov spectrum 



^k 



1/2 
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which can also be presented in the form 



-I 1/2 



2m 



2 4 

m c 



For the coefficient functions Uk and Vk, one has 



Ui. — 



2e, 



Vi. = 



2il 



which yields the properties 



2 2 1 



mc 

UkVk = - 7^ , 



2 I 2 



For the momentum distribution of excitations, one finds 



i/fe = < blbk > 



1 



e/3£fc - 1 ' 

while for the momentum distribution of atoms, one obtains 

rifc = < alttk > = (life + 't;fc)i/fc + vl 
For the anomalous average < a^aq >= < a^a^k >, one gets 

Note that the anomalous average is of the same order as n^. 
In the long- wave limit A; — > 0, one has Ek — chk, so that 



T^k ^ 



l3Ek 



mc 



1 



rik 



l3mc^ 1 
2 



1 /^gfc 

2 36 ' 



mc 



I3mc^ 
36 



2/3mc2 36 2 ' — - - 
Then the sum of the normal and anomalous averages gives 

, 1 (keT \ 

These properties are useful for calculating the momentum dispersion (68), for which one finds 

< p2 > = Y,{nkfnk - No lim{hkf < aktt-k > ■ (73) 
k ''^^ 



Introducing the notation for the mean kinetic energy 

rik 



^'=-"]v?^"'==2;;^y^"^(2^ 



(74) 



31 



and taking account of the limit 



lim(^A;) < aka^k >— —mkBT , 
fc— >o 



one comes to the expression of the normal density (66), which is 

2 E kill po , . 

Unfortunately, the mean kinetic energy (74) possesses an ultraviolet divergence, since at A; — > oo, 
one has Ek — {TikY /2m and nk — {mc/hk)^. Consequently, the normal density and, hence, the 
superfiuid density (65) are not defined in the Bogolubov approximation, while the condensate 
density (61) is perfectly defined. 

The message of this section is that a condensate and a superfiuid are not synonyms but are 
rather different things that should not be confused with each other. 



8 Collective Excitations 

Generally, collective excitations in dilute Bose gases can be of two types: elementary excitations 
and topological modes. Elementary excitations represent small oscillations above the ground state. 
They are described by the linear Bogolubov-de Gennes equations. In uniform gases, this yields the 
famous Bogolubov spectrum (see the previous Sec. 7). Elementary excitations leave the average 
spatial density unchanged. Contrary to this, topological modes are strong excitations essentially 
changing the spatial density distribution. These arc described by the nonlinear Gross-Pitaevskii 
equation. In uniform systems, the known example of a topological mode is a vortex. 

In nonuniform trapped gases, both these types of excitations are also present, though the 
situation becomes much richer. Elementary excitations in trapped gases can be subdivided in two 
kinds: continuous modes and discrete modes. 

When the excitation wavelengths are much shorter than the trap sizes, the excitation spectrum 
is a continuous function of its wave vector. This is why such excitations can be called contin- 
uous modes. The measurement of the k dependence of the excitation spectrum for the trapped 
condensate of ^^Rb atoms was accomphshed in the nice experiment by Steinhauer at al. [79]. 
The measured excitation spectrum was found to be in very good agreement with the Bogolubov 
spectrum in the local density approximation, even close to the long-wavelength limit of the region 
of applicability. The static structure factor S{k) was also measured and found to agree well with 
the Bogolubov theory. 

A Monte Carlo computation of the static structure factor S{k) of a dilute Bose- Einstein con- 
densate predicts a roton peak in S{k) occurring at k ~ Stt/o^, where Ug is the s-wave scattering 
length [169]. This peak is not steep enough to produce a minimum in the excitation spectrum e^. 

For excitations with wavelengths 27i/k which are comparable to the radius of the condensate, 
the excitation spectrum is characterized by discrete shape-dependent oscillatory modes. Because 
of this, such excitations can be termed discrete modes. For these modes, momentum is not 
anymore a good quantum number and consequently the elementary excitations are to be classified 
by the discrete quantum numbers related to the trap symmetry [2,4]. The excitation spectrum of 
discrete modes is mainly defined by the shape of the trapping potential, essentially not depending 
on whether bosons or fermions are trapped [170]. 
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Discrete modes have been observed in several experiments. Performing a small modulation of 
the trapping potential resulted in the observation of center-of-mass oscillations [171-174]. Higher 
multipolar-order surface oscillations were observed [175]. Scissors modes were also studied in 
experiments [176,177]. More details can be found in the review [2]. 

Theoretical calculations at zero temperature agree well with the observed frequencies of discrete 
modes [2-4]. Temperature dependence of the frequencies and of the corresponding attenuations, 
due to interactions between the condensate collective modes and noncondensed thermal atoms, 
have been theoretically analysed for uniform [178] and trapped [179-182] gases. Finite-temperature 
calculations are in good agreement with experiments for the lowest modes with axial angular 
momentum quantum numbers m — and 171 — 2 [183,184] and for scissors modes [185,186]. 

The standard way of treating Bose-condensed systems is based on the Bogolubov prescription 
ip = ipo + 'ip, with ipQ =< ip >, which breaks gauge symmetry. Then a Green function approach is 
involved for either real-time or imaginary-time Green functions, which are analytical continuations 
of each other on the complex-time plane [187]. Nonequilibrium processes can be considered by 
means of either Keldysh [91,188] or Kadanoff-Baym techniques [189,190]. 

Another approach is based on the Girardeau- Arnowitt canonical transformation 



[191,192]. This approach was followed in [181,193,194] and led to an accurate description of 
the lowest discrete modes [184]. Note that Bogolubov [70,71] also used a similar transformation 
ak — NQ^^'^alak, with the difference that here iVo is not an operator but the number of condensed 
atoms. In both methods, a well developed condensate is assumed to exist and in both the cases 
there appear the anomalous averages < a^ag >. The resulting equations are also similar in both 
these approaches. 

A very popular trick is when, after breaking gauge symmetry by means of the Bogolubov 
prescription, one omits all anomalous averages. One usually ascribes this trick to Popov. However, 
this is neither a correct approximation nor has any relation to Popov. Really, as is easy to infer 
from the original works by Popov [115,195,196], he never rejected the anomalous averages. Vise 
versa, he emphasized the necessity of accurately taking them into account in order to correctly 
calculate self-energies and to avoid nonphysical singularities after breaking gauge symmetry by 
the Bogolubov prescription. 

Moreover it is straightforwarded to show that, generally, the anomalous averages are of the 
same order as the normal ones, so that there is no any reasonable justification of omitting the 
former though leaving the latter [2,197,198]. For example, in the Bogolubov approach (see Sec. 
7) the normal average is 



where Sk is the Bogolubov spectrum and = {e^^'' — 1) ^, while the anomalous average is 




in which oq is the condensate operator and A^o is the condensate number-of-particle operator 
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These averages are closely related usually being of the same order. Thus, in the long-wave limit 
k ^ 0, one has < a^a^ >— mc^/Pel and < a^a^k >— —mc^/(3el.. That is, in this case, the 
absolute values of the normal and anomalous averages coincide at all temperatures below T^.. 
As soon as the anomalous averages appear, they must be considered together with the normal 
averages. There is no any justification of blindly throwing them out. 

Then one may ask: if the trick of throwing out anomalous averages is wrong, why, being used 
by so much people, it often yields reasonable results? For instance, it leads to accurate calculated 
frequencies of lower collective excitations in the temperature interval < T < 0.6Tc, though it is 
not correct for the m = 2 mode above 0.6Tc. The answer to this question is as follows. 

For a system with Bose-Einstein condensate, it is possible to employ a representation, based on 
mixed coherent states [60,199], where the gauge symmetry is not broken and, hence, the anomalous 
averages do not appear at all [125]. In this approach the field operator is presented as 

^(r,t)=r/(r,t)e^" + V^(r,t), (76) 

where 77 (r, t) is a coherent field and a e [0, 27r] is a random phase. The average of an operator A 
is defined by 

with p being a statistical operator. Because of the averaging over the random phase, all anomalous 
averages are identically zero, so that < ip >=< ip >= as well as < t/jip >—< ipip >— and so 
on. The equation for the coherent field r] = r]{r,t) reads 



drj \ 



V , (78) 



2m 

where U = U{r, t) is a trapping (external) potential, Ag = ATrh^ag/m, and p = p(r, t) is the density 
of noncondensed atoms. Equation (78) reminds a generalized Gross-Pitaevskii equation, but it 
should be stressed that Eq. (78) is an exact equation, with no approximation yet being involved. 
The density of noncondensed atoms can be found from the related single-particle Green function 

^(12) = -i < T^(l)^t(2) > , 

in which T is the chronological operator and each number j implies a set rj,tj. The evolution 
equation for G includes the corresponding two-particle Green function G2 and self-energy 



Sfl2) = A, 



25(12)|r/(l)|2 + i J (52(1113)^-^(32) d{3) 



This allows us to write the equation for G as 

Vl-U + p^ G(12) - J E(13)G(32) d(3) = 5(12) , (79) 

which is also an exact equation. These equations (78) and (79), being exact, at the same time 
contain no anomalous averages. More details can be found in Ref. [125]. Linearizing the coherent- 
field equation (78) with respect to the deviation 61] = r] — rjo, where rjo corresponds to the ground 
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state, it is possible to come to the Bogolubov-de Gennes type equations describing elementary 
collective excitations. 

In the stationary situation, when the external field U = U{r) does not depend on time, the 
coherent field can be presented as 



The eigenfunctions of Eq. (80) are termed topological coherent modes [200-203]. The mode, 
with the lowest energy Eq, describes the standard Bose- Einstein condensate, while the modes, 
with higher energies En > Eq, correspond to excited nonground-state condensates [200-203]. 
The generation of topological coherent modes can be realized by applying alternating fields with 
frequencies in resonance with the desired transition frequencies Umn = {Em — En)/h. The most 
known example of a topological mode is a vortex, which can also be excited by means of resonant 
external fields [204]. Reviews on vortices in trapped condensates can be found in [205,206]. A 
dipole topological mode was generated in experiment [207]. The properties of different topological 
modes were studied in several theoretical papers [200-203,207-221]. Such modes can also be 
generated by a resonant modulation of the scattering length [222-224]. The topological modes 
can be excited in single- well potentials, in optical lattices [225] , and in other multiwell potentials 
[226,227]. Adiabatic preparation and evolution of two tunnehng modes was investigated [228-230]. 
Thermodynamics of an ensemble of the coherent modes were described [231]. 

To stress more the principal difference between elementary excitations and topological modes, 
let us consider the single-particle density matrix pi, which can be presented as an expansion over 
its eigenfunctions \n >, called natural orbitals [65], so that pi = J2n^n\n >< n\, where the 
eigenvalue =< ajjOn > is the occupation number. In the case of a system with the usual Bose- 
Einstein condensate, the ground-state occupation number Nq r-j N becomes macroscopic, while all 
other levels, corresponding to elementary excitations, are not macroscopically occupied, so that 
Nn/N — > 0, as iV — > oo and n = 0. Contrary to this, for a system with several topological modes, 
several occupation numbers become macroscopic. This type of condensate is, sometimes, called 
fragmented [232]. Such a fragmentation in the spectrum of the single-particle density matrix 
should not be confused with a spatial fragmentation of a complex system [233]. Thus, a system 
with several topological modes is an example of the fragmented condensate in the sense of the 
spectral fragmentation [232]. It is also worth noting that the multimode fragmented condensate 
is, generally, nonequilibrium, since its realization requires an action of alternating external fields 
or a special initial preparation. 

9 Particle Correlations 

Trapped atoms usually form a very dilute gas, such that pa^ ~ 10^*^ — 10^'^. The quantity pa^ 
can be increased either by increasing the density p or the scattering length ag. The latter can 




Then, Eq. (78) reduces to the stationary eigenvalue problem 



(80) 



with the effective nonlinear Hamiltonian 



H[ip\ = - — V' + [/ + A, (a^oIv'I' + 2p) . 
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be varied by manipulating external magnetic fields close to Feshbach resonance [234]. These are 
collision resonances that occur when the energy of a colliding channel coincides with the energy of 
a long-lived bound state. When the colliding atoms and the bound state have different magnetic 
momenta, the resonance condition may be tuned via external magnetic fields exploiting the Zeeman 
effect [2,235]. If the parameter pa^ is essentially increased, then interparticle correlations could 
become important and even the whole picture based on the local interaction, involving just a 
scattering length, could turn invalid. Correlations become crucial, making the simple mean-field 
approximation unreliable, already at pa^ ~ 10^^ [54,236]. 

Feshbach resonances, first, were observed in Bose-Einstein condensates of ^^Rb, near the mag- 
netic field Bo = 160 G with width AS - -6 G [237], and of ^^Na, near Bo = 907 G with AB = 1 
G [238]. Then Feshbach resonances were realized in other bosonic atomic species [239-242], as well 
as in several fermionic atoms [243-246]. The behaviour of the scattering length near a Feshbach 
resonance at magnetic field Bo is given by the formula 



in which Ug = 05(00) is the off-resonance scattering length and AB characterizes the width of the 
resonance. Typical behaviour of the relative scattering length as{B)/as is illustrated in Fig. 1 for 
^^Rb, with Bo = 156 G and AB = -6 G, and in Fig. 2 for ^3 Na, with Bq = 907 G and AB = 1 
G. These figures follow directly from Eq. (81). In the case of ®^Rb, it was experimentally possible 
to change as{B) by about three orders of magnitude [237]. But for ^^Na, the variation of as{B) 
was limited by one order because of the strong enhancement of inelastic collisions close to the 
Feshbach resonance leading to an essential depletion of the condensate [238,247]. As is seen in 
Figures 1 and 2, it is possible not merely vary the scattering-length magnitude but even to change 
its sign. Such experiments were carried out with ®^Rb, when the sign of as{B) was fastly changed 
from positive to negative, which led to the explosion of the condensate [248,249]. 

When atomic interactions are so strong or density is so large that pa^ > 0.1, then it becomes 
crucial to take account of atomic correlations which are well characterized by the pair correla- 
tion function g{r,r') defined in Sec. 3. In many cases, it is sufficient to consider the uniform 
approximation g{r) depending on the difference of spatial variables. It may also happen that 
the contact potential is no longer representing interactions correctly and, instead of the contact 
potential proportional to 5(r), one should consider a more realistic form $(r). 

From the quantum scattering theory [250] it is known that, generally, the scattering ampli- 
tude is expressed through the integrals of the type / (/9*j(r)<I>(r)(/9„(r)(ir, where Lp{r) is a function 
characterizing the relative motion of two colliding particles and satisfying the Schrodinger equation 



in which rrieff = 17111712/ {mi + 777.2) is the reduced mass of the colliding particles with masses mi 
and 7772- For identical particles, with 7771 = 7772 = 777, the reduced mass 777e// — m/2. It is only 
for the low-energy Born approximation when the scattering amplitude contains just the integral 
/ $(r)(ir, which reduces to a scattering length. 

The simplest way of defining the pair correlation function would be by setting gir) = \ip{r)\'^, 
with ip{r) defined by Eq. (82), but with the boundary condition \(p{r)\ ^ 1 as r = |r| — > 00, 




(81) 



+ $(r) - E 99(r) = , 



(82) 



2777e// 
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instead of the standard normalization condition. Such a definition of g{r) was suggested by Bo- 
golubov [70] and employed for describing quantum crystals [251-257], quantum liquids [258,259], 
melting-crystallization phase transformation [260], and, generally, for considering quantum sta- 
tistical systems with strongly singular interaction potentials [261]. The parameter E in Eq. (82) 
has to be chosen from some additional physical conditions. Thus, it is natural to assume that the 
pair correlation function g[r) = |<^(r)p possesses a maximum at Tq where the interaction potential 
$(r) has a minimum. Then one should set E — $(ro). For a short-range repulsive potential, one 
may set E = 0. 

For dense systems, it is possible to develop more refined ways of constructing the pair corre- 
lation function g{r). It is also possible to find g{r) from experiments [62,262] by measuring the 
structural factor. 

Including interparticle correlations from the very beginning implies that the interaction po- 
tential $(r) should enter everywhere together with the pair correlation function g{r), forming 
the smoothed potential $(r) = ^{r)g{r), which plays the role of a pseudopotential. The latter 
is always finite and integrable, since, when $(r) becomes singular, the correlation function g{r) 
tends to zero, so that the smoothed potential is always nonsingular, which removes the ultraviolet 
divergencies appearing in the standard consideration starting with a simple mean-field approxi- 
mation, such as the Hartree-Fock approximation. An approach, involving at all iteration steps 
only the smoothed potential $(r), was, first, suggested in [263], which has been developed into a 
refined correlated iteration theory [60,264,265]. Here we present only the final set of equations of 
this theory. 

For the single-particle Green function 

^(12) = -i < 7^(1)^(2) > , 

we may write the Dyson equation 

(7(12) = ^0(12) + J (70(13) [E(32) - Eo(32)] d{3) , (83) 

in which 6*0(12) is any trial Green function associated with the self-energy Eo(12). Here the same 
abbreviated notation as in Sec. 8 is used. The self-energy 

E(12) = ±« / ^(13)D(132) d(3) (84) 

is expressed through the smoothed potential 

^12) = ^(12)$(12) , (85) 

where 5'(12) is the pair correlation function, and the doubling function D(123) is given by the 
form 

D(123) = YDo{123) , (86) 
where the starting doubling function is 

Do(123) = (5(13)G(22) ± (7(12)5(23) , 

and the operator Y is defined below. For generality, Bose as well as Fermi statistics are implied 
here, with the upper sign being related to Bose, while the lower sign, to Fermi statistics. 
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An iterative solution of Eqs. (84) and (86) follows the scheme 

Dk Sfc+i — >• Ifc+i — > -Dfe+i , (87) 
where the /c-order approximation for the operator Y is 

k 

+ (88) 

and the action of the operator on any three-point function /(123) is defined as 

yfc/(123) = [l-5(12)]/(123)+ 

+ / G(14) [^(52)/(527)G'(76) - G(22)G(56)5(67)] d(4567) . (89) 

At any approximation order, particle interactions enter always in the form of the smoothed po- 
tential (85) renormalized by the pair correlation function 5'(12). Thus, particle correlations are 
taken into account from the very first step of the correlated iteration theory [60,264,265]. 

An analogous procedure can be followed for finding the response functions [60] . For instance, 
let us consider the response function 

where 5/i(l) is a variation of the chemical potential. The poles in the energy representation of the 
function x(12) = x(112) describe the spectrum of collective excitations. The response function 
(90) can be presented as the action of an operator Z on the initial function 

Xo(123) = ±zG(13)G(23) , 

corresponding to the random-phase approximation, so that 

X(123) = Zxo(123) . (91) 

An iterative solution of the latter equation is realized according to the scheme 

Sfc ^ Zfe ^ Xfc , (92) 

in which the A;-order approximation of the operator Z is 

k 

ife = l+^f^ (93) 

with the action of prescribed by the rule 

4/(123) ^ / G(14)G(52) /(673) d(4567) . (94) 
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Here again, we never meet the bare interaction potential, since the self-energy Efe(12) contains 
always only the smoothed potential (85). 

Wishing to take an explicit account of the arising Bose-Einstcin condensate, we may proceed as 
in Sec. 8, shifting the total field operator as in Eq. (76) and separating an effective wave equation 
for the coherent field, as in Eq. (78), from the equation for the Green functions of noncondensed 
particles. However, a more general way is to present the total Green function as the sum 

G(12) = -^r/(l)r/*(2) + G'(12) . (95) 

Note that the shift (76) is sufficient, though not necessary, for the representation (95). A system, 
characterized by the Green function (95), contains a coherent part, related to the coherent field 1], 
and a noncoherent part of noncondensed particles corresponding to G. The existence of a coherent 
part can be detected experimentally [266]. 

The full many-body theory, accurately taking into account particle interactions and correla- 
tions, will become necessary as soon as the parameter pa^ is not anymore small. As is explained 
at the beginning of this section, this parameter can be made sufficiently large, up to pa^ ~ 1, by 
tuning atomic interactions with an external magnetic field close to a Feshbach resonance. In the 
vicinity of the latter, there exist bound states, when atoms can form molecules, whose interactions 
are long-ranged, also requiring to accurately take into account particle correlations [267] . 

10 Molecular Condensation 

Since there exist bound or quasibound states for colliding cold atoms, it has been an impor- 
tant challenge to produce ultracold molecules that could form Bose-Einstein condensates. Such 
molecules would behave as bosons, independently of whether they are composed of bosons or 
fermions. 

One route for producing ultracold molecules is to form them from Bose atoms in the condensed 
state. Thus, a two-photon stimulated Raman transition in the condensed ^^Rb was used to 
produce ^''Rb2 molecules in a rotational- vibrational state [268], and ultracold molecules of ^'^Na2 
were formed [269] through photoassociation of the condensed ^^Na. The molecules of ^^Na2 were 
also produced by means of the Feshbach resonance technique [270], which was used as well for 
producing the molecules ^^^Cs2 [271]. Employing the Feshbach resonance was found to be the most 
convenient way of creating cold molecules, and this technique was used in all experimental works 
cited below. Molecules are formed by sweeping an external magnetic field through the Feshbach 
resonance, adiabatically converting atoms to molecules. 

The detection of whether the created molecules form a Bose-Einstein condensate is a demanding 
problem requiring to probe coherence properties of the assembly of atoms plus molecules. It was 
found that the molecules ®^Rb2 do form the condensate [16] as well as the molecules ^''Rb2 [272]. 
The molecules made of cold Bose atoms are in highly vibrationally excited states and usually 
undergo fast decay during lO^'' s. 

Molecules can also be composed of cold Fermi atoms by tuning the atomic interactions with a 
varying magnetic field close to a Feshbach resonance. It is possible to vary the s-wave scattering 
length tts from positive to negative over many orders of magnitude. This opens up the possibility 
of investigating the crossover from a molecular condensate to atomic superconductivity [273,274]. 
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The condensate limit occurs for a positive scattering length > 0, much smaller than the inter- 
atomic spacing, so that pa^ <^ 1. Then the fcrmions form weakly bound dimers of a size about a^. 
The superconducting limit is expected to occur for a negative scattering length < 0, when the 
distance between fermions in a superconducting pair is much larger than the interatomic spacing. 
Then the Fermi momentum A;^ is such that A;f|Os| <S 1. The Fermi temperature of a harmonically 
trapixxi gas is Tp ~ (/ia;o/A;B)(3iV)^/'^. The superconducting transition is expected to happen 
when T/Tp ~ 10"^ - 10"^ Thus, if Tp ~ lO^^ K, one needs to reach the temperatures T ~ 10"^ 
K, which is achievable in modern experiments. 

The molecules ®Li2, composed of cold Fermi atoms ®Li, where formed near a Feshbach resonance 
[275-277], being rather long-hved, with a lifetime around 1 s. The molecules ^°K2 hve shorter, 
about 10~^ s [278]. It has been possible to convert to molecules up to 80% of atoms. 

Molecules, formed from Fermi atoms, behave themselves as bosons and, hence, can be con- 
densed at sufficiently low temperatures. Thus, the molecules ^Li2 condense at = 6 x 10~^ K; 
below this temperature condensate fraction of up to 90% has been observed [17,279-282]. The 
molecules ^'^K2 can also be condensed, though their condensate fraction being about 10% [283,284]. 
It is especially interesting to study the crossover region around a Feshbach resonance, going from 
positive scattering lengths to negative ones. Then the molecular condensate transforms into a 
degenerate Fermi gas of atoms. This crossover region was investigated for [280-282] as well 
as for '^°K2 [284]. In the latter paper, condensed molecules of K2 were observed in the region 
where the scattering length is negative and no stable ^°K2 molecules would exist in vacuum. It 
was claimed [284] that those condensed pairs were not bound into molecules, but merely moved 
together in a correlated fashion, similar to Cooper pairs of electrons in a superconductor, hence 
the claim was that the superconductivity of '^^K atoms had been observed, although the exact 
nature of these pairs remained unclear. However, theoretical calculations showed [285] that these 
pairs are not the Cooper pairs but are the same tightly-bound molecules as those existing in the 
condensate region, where > 0. The existence of such molecules in the degenerate Fermi-gas 
region, where < 0, is due to the interaction of these molecules with the unbound atoms of 
the degenerate Fermi gas. This conclusion was con&med by the experiment [282] with condensed 
molecules ®Li2, where the large condensate fraction, up to 80%, was observed in the degenerate-gas 
region, where no stable ^Li2 molecules could exist in vacuum. The existence of these molecules 
has been interpreted as been stabilized by the presence of the Fermi gas of unbound atoms [282]. 

The situation, when bound states, or bound clusters, that cannot exist in vacuum, nevertheless, 
can perfectly survive inside a sea of their unbound constituents is not unique for trapped atoms 
[122]. The same picture happens for hadrons that are stable below the point of deconfimenent 
but which would not be stable above this point, if they would be placed in vacuum. Nonetheless, 
a fraction of hadrons can perfectly survive above the deconfinement point provided they are 
immersed into quark-gluon plasma [121, 122]. 

Concluding this review, I would like to repeat what was said in the Introduction. The physics of 
trapped cold atoms is a very wide and fast developing field. A number of interesting and important 
subjects are not touched here, such as spinor condensates, dipolar condensates, multicomponent 
condensates, and condensates in optical lattices. I have concentrated in the present review on the 
most fundamental, to my understanding, problems. And I hope to cover other topics, not touched 
here, in another publication. 
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Figure Captions 



Fig. 1. Relative scattering length as{B)/as for ^^Rb close to the Feshbach resonance at 
Ba = 156 G with AB = -6 G. 



Fig. 2. Relative scattering length as{B)/as for ^^Na close to the Feshbach resonance at 
Bq = 907 G with = 1 G. 
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